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Abstract: We calculate the high energy behavior of quark-antiquark exchange in 7* 7* 
elastic scattering by summing, to all orders, the leading double logarithmic contributions of 
the QCD ladder diagrams. Motivation comes from the LEP data for cr]^^ which indicate 
the need for secondary reggeon exchange. We show that, for large photon virtualities, this 
exchange is calculable in pQCD. This applies, in particular, to parts of the LEP kinematic 
region. 
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1. Introduction 

The 7* 7* collision at high energies is the unique laboratory for testing asymptotic prop- 
erties of perturbative QCD. The virtuality of the photons justifies the use of perturbative 
QCD, and modern electron positron colliders (LEPII, a future linear collider NLC) allow to 
measure the total cross section of 7* 7* scattering at energies where asymptotic predictions 
of perturbative QCD can be expected to set in. The dominant contribution to the process 
is given by the BFKL Pomeron which gives rise to a cross section strongly rising with 
energy crj^^i ~ s"^^''-'. Here ap(0) is the Pomeron intercept which, in leading order and 
for realistic values of the photon virtualities, lies in the region ap(0) ~ 0.3 — 0.5. 

There is little doubt that the pomeron will dominate at very high energies, and it is 
expected to be a main contribution at any future linear collider. At present, however, the 
only source for experimental data on photon photon collisions is LEP ^. These data are 
at energies which cannot be considered as asymptotically large, and it has become clear 
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that at LEP energies, the cross section is not yet dominated by the pomeron |Q, |5|, |6|, 
1^, ^, 1^ . The data rather indicate the necessity to include, in the theoretical description, 
several corrections. Perturbative corrections are due to the quark box (often referred to 
as QPM contribution); recently first 0{as) corrections to the quark box have been 
considered. NLO corrections to the BFKL Pomeron are not yet fully available: whereas 
the NLO corrections to the BFKL kernel have been completed [0], those of the photon 
impact factor are not yet finished |12, Nonperturbative contributions include the soft 
Pomeron (in the low-Q^ region) and the exchange of secondary reggeons: the exchange of 
/o (flavor singlet) or Aq, A2 (flavor nonsing let) In this paper we address the latter 

contributions, the corrections due to secondary reggeons. 

In hadron hadron scattering, secondary reggeons denote the exchange of mesons and 
are of nonperturbative nature. This may change, however, if we replace one of the hadrons 
(or both) by a virtual photon. We know from the 7*p process in deep inelastic electron 
proton scattering that the rise of the total cross section with energy becomes substantially 
steeper than in proton proton scattering; this observation has lead to the notion of a 'hard 
Pomeron'. What makes this hard Pomeron particularly attractive is that its energy depen- 
dence becomes calculable in pQCD. In 7*7* scattering, the high energy behavior should be 
described by the BFKL Pomeron, i.e. pQCD provides an absolute prediction for this pro- 
cess. In analogy with this, one might expect that also the secondary exchange may become 
accessible to a perturbative analysis, if we replace one or both external hadrons by virtual 
photons. If so, meson exchange will be modelled by the exchange of qq ladders and its 
prediction for the energy dependence may be tested in the corrections to BFKL exchange 
in 7*7* scattering. It is the purpose of this paper (and a forthcoming one) to investigate 
this possibility more closely. As a start, we present a perturbative QCD calculation of 
the even-signature flavor nonsinglet contribution to the total 7* 7* cross section. This is 
done by summing over gluon ladders with two quarks in the t-channel. The flavor singlet 
case appears to be more complicated since it involves the mixing between quark-antiquark 
and two-gluon states in the t-channel (where the two-gluon state has a helicity content 
different from the so-called nonsense- helicity configuration in the BFKL Pomeron). The 
flavor singlet exchange is expected to be somewhat larger than the nonsinglet one |jl5|, and 
we will come back to it in a forthcoming paper. 

One of the striking differences between gluon exchange in the BFKL calculations and 
quark-antiquark exchange is the appearance of double logarithms |16, 17 1. As result of this, 
the intercept of the qq-system is of the order Wq'' = ^/ const Os (as opposed to oJq^^^ = 
const as in the single logarithmic high energy behavior of the BFKL Pomeron), and its 
numerical value can be expected to be large. In fact, for qq scattering it is known |14| 
that the cross section goes as ~ s'^°~^ with wo = ^2as C^/vr ~ 0.5. It is remarkable 
that this intercept obtained in pQCD is very close to the nonperturbative one known from 
phenomenology. For 7* 7* scattering we obtain the same result for the intercept. 

Another crucial feature of the double logarithmic calculation is its dependence on the 
infrared region. In principle, the ladder graphs to be summed are infrared safe. However, 
quite in analogy with the BFKL approximation, the contribution from small momenta is 
not believable, and one has to introduce a cutoff scale, /io < IGeV, which separates the 
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infrared from the large momentum region. Taking the virtuahties of the external photons 
to be sufficiently large and then considering the high energy limit one observes that, for 
not too large energies, the result is independent of this infrared cutoff: a natural infrared 
cutoff for the transverse momenta inside the ladder appears, which is of the order Q'^/s 
and decreases with energy. Thus with increasing energy this cutoff eventually reaches the 
scale yUgj and the integration starts to penetrate the nonperturbative domain. In this region 
perturbative QCD cannot be trusted any more. Therefore, for energies s larger than Q^/ijCq 
one has to limit the transverse momenta to be larger than //q, and the perturbative high 
energy behavior of the qq exchange starts to depend upon the cutoff. This dependence 
turns out to be rather weak: it mainly resides in the pre-exponent and not in the exponent 
of the energy. In our calculation we will study this dependence in detail: we will compute 
the sum of double logarithms 



n>l 



^a„arMln2^) . (1.1) 



When ^ reaches the cutoff /Xq another large logarithm appears, InQ^/^Ug, which we will 
include into our analysis: 



n>l ^ 



a„o In^" 7^ + a„i In^""! ^ In ^ + ... + a„„ In^" 9l ) . (1.2 



In particular we will study the role of /Xq in the determination of the pre-exponential 
behavior of the high energy asymptotics. 

The task of the double log resummation in 7* 7* collision can be attacked by several 
methods. First, we follow the original paper Ref. jl^ in which the question of double log 
resummation e^e~ annihilation was addressed first. By direct summation of the ladder 
graphs this method leads to a Bethe-Salpeter type equation for the amplitude. Compared 



to the original work of Ref. |16] the case of 7*7* is slightly more complicated, since the 
additional scale is involved. This complication results in a more sophisticated structure 
of the solution. We apply this method to ladder diagrams which are the relevant diagrams 
for the even-signature exchange. One of our main results is the discovery of the 'hard 
region' /Xg < Q'^/s in which pQCD is reliable. 

The second method uses the infrared evolution equation (IREE) for the amplitude 



in the Mellin space. IREE was first derived in Ref. |17] in application to quark-quark 
scattering, and, more recently, has been applied to the calculation of the small-x behavior 
of flavor nonsinglet structure functions in Ref. and to polarized structure functions 
in Refs. [|^, 19 1. Within the double log accuracy IREE traces the dependence of the 
amplitude on the infrared cutoff fi of the quark transverse momentum in the ladder. The 
scale fj, is auxiliary in this method. We will see, however, that when identified with the real 
scale of nonperturbative physics, ^q, it leads to the same answer as the linear equation of 
the first approach. Use of IREE allows to discuss not only the even signature case (with 
ladder diagrams), but also the odd signature case (with non-ladder diagrams). 



A third way of handling the quark-antiquark exchange has been described in [EG, 21 



using the partial wave formalism and the notion of the reggeon Green's function, a linear 
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equation a la BFKL has been found for the fermion-antifermion scattering amphtude. 
In 1 21 1 the kernel has been shown to have the same conformal properties as the BFKL 
kernel. However, the application of this formulation to our problem is not completely 
straightforward: the coupling of the quark-antiquark ladder to the external photons involves 
a careful handling of the double logarithmic integrations. We will show how our result for 
the hard region can be obtained from the reggeon Greens's function. 

The structure of the paper is following. In the next Section (2) notations are intro- 
duced, and the quark box diagram is discussed. In Section 3 we sum the gluon ladder and 
derive a Bethe-Salpeter type equation for the amplitude. This equation is solved in Section 
4. We also discuss the solution and present the high energy asymptotics. In Section 5 we 
rederive the same result using the method of the infrared evolution equation (IREE) . The 
approach based on the reggeon Green function is discussed in Section 6. Some numerical 
estimates are given in Section 7. Section 8 contains our conclusions. Two appendixes 
contain some details of our calculations, supplementing sections 4 and 5. 



2. The quark box 

We begin with the lowest order diagrams for the scattering amplitude T'^ of the elastic 
7*7* scattering process. We restrict ourselves to the forward direction i = 0, and for 
simplicity we first take the virtualities of all external photons to be equal. 

we restrict ourselves to the high 



The exact computation can be found in Ref.|22 



energy behavior. The lowest order consists of the three fermion-loop diagrams (Fig. |T], a - 
c); at high energies only the planar box (a) and its s <-> u cross partner (b) contribute to 
the leading double-logarithmic behavior. 
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Figure 1: The Born level diagrams. 

The calculations below are done using the Feynman gauge. Our method of extracting 
the double-logs will be close to the original paper We introduce the notation (Fig. |l]) 

pa = -Q^: pI = -Q^; {va+Pb? = s- x = Q^/s. . (2.1) 

In our Sudakov decomposition k = Pq — ap + k± (with k'j_ = —}?) the light cone vectors 
p, q are defined through 

PA = P - xq; pB = q - xp; p^ = q^ = 2{pq) = s. (2.2) 
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Transverse polarization vectors are defined by e^^. = (0, l,ibi,0), the longitudinal 
ones by eL{pA) = {p + xq) = {pA + 2xq), eLips) = {p + xq) = 

{pB + 2xq). 

We first analyze the quantum number structure of the t-channel. We Fierz-transform 
the upper part of the trace expression: 

Y{^+i>A)l^ = {{k+pATg^' + {k + pArg"' - {k + PAYg^nip + ie'"'^''ik+PA)al5lp. 

(2.3) 

As discussed in [^], the vector current generates the even signature A2 exchange, the 
pseudovector current the Ai exchange. The absence of the scalar, pseudoscalar, and tensor 
currents seems to indicate that, in 7*7* scattering, there is no scalar or pseudoscalar 
(vr-like) exchange; note, however, that the first gluon rung inside the quark loop might 
contain the axial anomaly, i.e. the coupling of the t-channel pion-like state to two gammas 
is nonzero. In our subsequent analysis we will restrict ourselves to the vector contribution. 
The pseudovector exchange as well as the pseudoscalar exchange decouple from the helicity 
conserving scattering, i.e. is of no interest for the total cross section. 

Let us now turn to the quark loop. In order to find a double-logarithmic contribution 
we need to find, in the trace expression in the numerator, terms proportional the leading 
power of s and to k'j_. For the numerator of the planar box diagram with transverse 
polarization we obtain: 

tr {A)^m{^-^B)^ {B') HiA') {i>A + ^PA-PBtr {A) (B') , 

(2.4) 

where, on the rhs, we have made use of the fact that, in the high energy limit, the leading 
power of s is due to the product pA ■ Pb- For the helicity conserving scattering we obtain: 

tr{^iA)^mi^-^B)^{B')^^{A'){^A + f)) ^2skl x 

(2.5) 

[e{A) . e{A') e{B) • e{B') + e{A) • e{B') e{B) • e{J^) - e{A) • e{B) e{^) ■ e{B')\ . 

For the helicity conserving case (which we will need for the total cross section) the last two 
terms cancel. In the following, therefore, we shall restrict ourselves to the first term of the 
rhs of (12.51). For later convenience we define the trace factor: 



TTT = 4 AT, al^ Fns e{A) ■ e{A') e{B) ■ e{B'), (2.6) 

where 

quarks ■' quarks 

denotes the projection on the fiavor nonsinglet t-channel (for the flavor group SU{Nf)), 
and Cq stands for the electric charge of the quark, measured in units of e. 

Diagram Fig. ||b will be obtained by simply substituting, at the end of our calculation, 
s ^ u. The nonplanar diagram Fig. ||c has no contribution proportional to sk'j_ and will 
be neglected. 
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As to the integration over the momentum k, it will be instructive to first follow |16|: 
taking into account the k'j_ factor from the numerator, we first do the transverse integration 
(formally by replacing one of the exchange propagators by a 5 function), and are then left 
with the a and /3-integrals of the type 



da f dp 

^ J T 



(2. 



which are restricted to lie in the region: 

X < a < 1, X < P < 1, fil<saP = P<s (2.9) 

where fio denotes the momentum scale which separates the infrared (nonperturbative) 
region from the hard region. In the following, however, we will use, as integration variables, 
P and k"^, and the limits have to be derived from ( ^.91) . We have to distinguish between 
the two kinematic regions which we denote by '+' and '— ': 

/+ : < ^ (2.10) 

I- : ^ < /xg. (2.11) 
In the first region, /"'", we have the following limits of integration 

' ^ f'!^ = A. (2.12) 
In the second region, , we have two separate contributions: 

^o/Q' dp r^^ d^ dp j-f^' d^ _ 1 / ^2 _±_ _ i^^i^9l 

(2.13) 

One easily verifies that at = ^ the two results coincide. All other photon po- 
larizations, at high energies, are suppressed compared to the helicity conserving case of 
transverse polarization. For example, the case TL TL goes as const, the case LL — > LL 
as 1/s. In the following we will restrict ourselves to the leading polarization, to TT TT. 
Combining (U) and (|2l|), (|2l3|) we write our result for the planar box diagram in the 
form 

The nonplanar box in Fig. |l|b is obtained by substituting s — > u, and in the sum of 
the two diagrams the obtained results stand for the even signature A2 exchange. The total 
cross section for 7* 7* (averaged over the incoming transverse helicities) follows from 

a^*^* = ^ImT c.-^ (2.15) 
s s olns 
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(with the last approximate equahty being vahd in the high energy approximation only) 
and has the form 

i 2 In if Mn < — 1 



■ Mo 



It is not difficult to generalize our analysis to the case of unequal photon masses, 
and Instead of the scaling variable x = Q'^/s we now have the two variables Xi = Qf/s 
{i = 1,2). The integration region ( p.9| ) is replace by 

X2 < a < I, xi < (3 <l, ^Q<saj3 = }?<s. (2.17) 

The two kinematic regions ( p.lOD , ( 2.11| ) become 



/+ : 9jM (2.18) 

r: 9M< (2.19) 
and the result ( |2.14| ) for the planar box diagram takes the form: 



3. Linear equations for the ladder 

We now turn to higher order corrections to the quark loop diagrams. 

For the gluons we will use Feynman gauge, and we make use of the discussion given 





m 



16, 17, nS, 19, B3]. When selecting those Feynman diagrams which in the high 



energy limit contribute to the double logarithmic behavior, one has to distinguish between 



even and odd signature. For the case of quark-quark scattering it has been shown [17| that 
the even signature amplitude is described by ladder diagrams, whereas the odd signature 
amplitude has also non-ladder graphs. In this and in the following section we derive and 
discuss a linear integral equation for the even signature amplitude; so we can restrict 
ourselves to QCD ladder diagrams. The odd signature case can be derived from the IREE 
(to be discussed in section 5). 

Let us consider the diagram with one gluon rung (Fig. ^). We begin with the trace 
and look for terms which result, for each of the two loop momenta ki and k2, in a factor 
kf. Writing down the trace we find, for the lower cell, the product ...'^^ ^2 ^{B) (^ — 
^b) HB') h ■■■ which can be rewritten as 2A;|_l ^2 H^) (^ - ^b) H^') h + ■■■ where 
the remainder (indicated by dots), for the helicity conserving case e{B) = e{B'), does not 
contribute to the double logarithmic approximation. The result for the trace therefore 
becomes 

Askj^kl^eiA) . e{A') e{B) • e{B'). (3.1) 
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\ 

(a) (b) 
Figure 2: The first gluon rung (a). The ladder (b). 
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The k'^ factors will cancel against exchange propagators; together with color and vr factors 
we obtain the factor 

A = (3.2) 

for the gluon rung. 

Turning to the momentum integrations, we, once more, first follow the procedure 



outlined in [16| and integrate over the transverse momenta. The remaining a and /? 
integrals are of the logarithmic type (|2.5| ), and the variables are ordered according to 
X < 02 < cti < 1, X < /?! < /32 < 1, and their products are restricted by < 
saiPi = kf < s. In what follows, however, we choose the variables /3j, kf = kf; the 
integrals to be done are 

dPi r dkf r d(32 

~K J ^ J ^ 

with the integration region 



dkl 

^2 



(3.3) 



X 



— < /3i < /32 < 1, max{j32Q'^, //q) < ^2 < P'^s, max{^ — , /Xg) < kf < f3i s . 
s P2 

(3.4) 

As in the one loop case, we shall see that the two possibilities ( 2.10| ) and ( 2.1l| ) will have 
to be distinguished. We start with the integration of the upper cell, keeping the variables 
of the lower cell, /?2, ^21 fixed. The restrictions (3^) imply two distinct regions (Fig. ^): 



2 , ^2 



^0 < 



S/32 



k ? 



In the first case we have 



/32 



dpi 

2/s /^l Jl3ik^/I32 



, s/?2 , Sf32 

in — TT in 



(3.5) 
(3.6) 

(3.7) 
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whereas the second one leads to a sum of two contributions: 



Pi 



Mo 



In 



^2 



+ 



/32 



dkl 



dpi 



In' 



2^! 



^ 1 ^2% /?2S 



(3.8) 



For the remaining integrations over (52 and k"^ (Fig. ^) we not only have to observe the 
distinction ( |3.5D , ( |3.6[ ), but also to return to the two cases defined in ( 2.1C| ), ( |2.11 ). In the 




(a) ^^l<^ 



(b) ,l>^ 



Figure 4: The integration domains for one rung: lower loop. 

first case (i) we find (Fig. ^) that the line /cl = P2 always lies above the line P2 /^o ^/Q^ 
and max{P2Q'^, fJ-o) = /?2Q^- Therefore in the whole integration region (shaded area in 
Fig. ^) we have A;|/s > P2fJ-'o/Q'^, and the first condition ( |3.5D is always fulfilled. This 
leads to 

d(32 f^'' dkl B2S B2S I . s 
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For the second case (ii) in ( |3.6D the situation is shghtly more compHcated (Fig. ^p): 
the region of integration (shaded area) consists of three pieces (denoted by A, B, and C). 
In C the first condition ( |3.5| ) holds, whereas in B and C the second one (^^) is fulfilled. 
For region A the result is: 





dP2 




dkl / 




132 







1 ln4 4^ - i ln2 4- ln2^ + - In 4, In^ ^ - - In^ ^ . (3.10) 
12 Q2 2 Q2 ^2 ^ 3 g2 ^2 4 ^2 V 7 



For region B we obtain: 



f^l/Q^ /?2 7^32 Q2 k1 \2 \ kj ^il) III kl 



= i in3-f- In^ + - ln2-f- In^^ - In In^ ^ + - In^ ^ . (3.11) 

3 g2 ^2 4 Q2 ^2 Q2 ^2 12 > 

Finally, region C: 



P/s h y/32SMfVQ2 ^2 <3 ^2 4 Q2 



Adding up these three contributions we find 

712 + 2 — + o -2 In — - In -2 In^ — 



At the point = ^ both results (|3.9| ) and ( ^.13 ) coincide. 

After the analysis of the one rung ladder diagram it is not difficult to see the general 
pattern and to derive an integral equation. The general diagram with n rungs is illustrated 
in Fig. |2|b. ^From the trace expression in the numerator we obtain, for each rung, a factor 
kf^ A, where the momentum factor cancels against one of the exchange propagators of the 
cell labeled by "i". Following the strategy adopted for the one loop diagram and for the 
two loop case, we arrive at the a and /3 integrals of the logarithmic type ( |2.8[ ), with the 
ordering conditions 

X < an < a„_i...Q!i < 1, X < j3i < (32 < ■■■ < (3n-l < /3n < 1 (3-14) 

and the additional constraints 

fiQ < s ai Pi = kf < s . (3.15) 

Translating to the variables f3i and kf = kf, the limits of integration of the kf become, in 
analogy with (^3): 

Pn-i k'^ 

max{PnQ'^, /^o) < kn < PnS, max{-^ — -, /Uq) < A;^„i < Pn-i s 

max{^^, fx^) < kl < Pis. (3.16) 

P2 
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Starting from cell 1 at the top, we now work our way down along the ladder. We first do 
the integrals of cell 1, keeping fixed the variables of all the other cells. In doing so, we are 
led to distinguish the two cases ( ^ ) and ( |3.6| ); the results are given in ( p.7D and (^), 
which we denote by Af{l32, and ^^(/32, A;!), resp. The following integration of cell 2 
is done as described before: the only difference is that (2.10) and ( 2.11| ) are replaced by 

/3+: /.g<^; (3.17) 

I^: ^<f^l (3.18) 

As before, the '+' region receives contributions only from the '+' region (FigJ^), whereas 
the ' — ' region has contributions from both region A and B and from region C (Fig.^b). 
Denoting the results by A^(/33, k^) and (/?3, A;!), we have, symbolically, 

A+ = K++ A+ (3.19) 

and 

A- = ® A'l + K— (g) A^ . (3.20) 

Note, in particular, that the evolution of A'^ decouples from A^ . The same pattern now 
repeats itself for the integration of cell 3 etc, and we can write down a recursion relation. 
Instead, we define (/?, k'^) = A^ (where = 1) which satisfy the coupled integral 
equations (Fig. P): 



and 



^+(/3, e) = l + X ^ k'^) (3.21) 



A-{P, k') = 1 + X 1' ^ r ^A^iP', k'^) + 

(3.22) 



A/ / -j^A~{P',k') + X A~(p',k'). 

Jq2/s p J til Jpi4/k^ P Jp'k-i/p K 



A(k,p) 



k,|3 



A(k;p') 




Figure 5: The Bethe - Salpeter equation. 
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The first equation applies to the region , the second one to 7^ . They are defined by 



A ■ 



1-^0 < 



k'^ 2 
< ^0- 



(3.23) 



(3.24) 



Finally, the amplitudes for the photon-photon scattering are obtained from A by 

± 





subtracting the terms = 1 and by putting k"^ = and P = 1: 



tHq',s) = [a^{i, Q^)-i]^. 



(3.25) 



Note that in this limit the regions coincide with and in (2.10) and (2.11). 

The full scattering amplitude is obtained by adding the twisted (with respect to s u 
crossing) fermion loop. 

Before we finish this section we again briefly present the generalization to the case of 
unequal photon masses. With the main differences being the lower limits of the a and P 
variables (analogous to ( 2.17| )), the analysis for the one gluon rung and for the derivation 
of the integral equations can be done in exactly the same way as for the equal mass case; 
the final results for the integral equations are obtained from ( 3.21 ) and ( 3.22| ) by simply 
substituting Q"^ ^ Q\: 



A+{p,k^ 



1 + A ^ 



^ dp' f^'' dk' 



Ql/s Jl3'k^//3 



k'2 



'T^A^{P',k'^) 



(3.26) 



and 



A-ip,e) = i + x f 



%A^{P', P) + 



A 



j/s P' Jp'^.l/Ql k"^ 



(3.27) 



/3 dp' r^'f^o/Qldk' 



^■nrA-iP', k'^). 



k 



The same replacement applies to the kinematic regions (1^), (|3^ ). At the end we put 
A;2 = Ql and P = 1: 



T^iQlQls) = [A^{l,Ql) - 1] 



(3.28) 



For the remainder of our paper we will stay with this general case of unequal photon masses. 



4. Solution of the linear equations 

The structure of the two equations, ( |3.26|) and ( f3.27 ) defines our strategy: we first solve 
the equation for A"*", Eq. ( |3.26| ), and we then use the solution as an inhomogeneous term 
in the equation for A~ , Eq. ( |3.27D . 
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4.1 Solution in the region 

In the region where Eq. ( p.21 ) holds we define the new variables: 



e = ln(/3sA2); 7? = ln(/?s/Q2). 

In these new variables the equation ( 3.21| ) can be rewritten 



A+iC, 77) = 1 + A r dr? dWiL fj) . (4.1) 
Jo Jo 

When differentiated twice Eq. (^) reduces to 

A solution to Eq. (|4.2| ) can be found in the following form 

A+{(,V) = I ^Mz)e^^ + ^^'\ (4.3) 
Jo+ 2 7rz 

with the integration path closed around the point z = Q. 

The function (t>+{z) is uniquely determined when (|4.3D is substituted into Eq. (0). 
The result is 

ct>+{z) = (4.4) 
z 

and 



A+{^,1^) = iJ^/IX^] . (4.5) 



4.2 Solution in the region 

In the region we have the integral equation ( |3.22| ) which couples the functions A'^ and 
A-. 

Again we define a new variable: 

i' = i- Lo = HQl/lA) ■ 

In the variables , rj) the equation ( 3.22| ) can be rewritten 



v) = 1 + ^ r df) f A+(f , 7?) 
Jo J-Lo 



+ A / dn r d^'A~{C', f)) + \ r df] I di'A-{i', fj) . (4.6) 
Jo Jo Je Jo 



2 



Note that, at f = 0, ^"(O,??) coincides with ^+(Lo,?7): = f = i - Lq = In 

denotes the point where the two regions l\ and touch each other. 
After double differentiation ( [4.6D reduces to 
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For the solution we use the same ansatz as before: 

dz . , , 



0^ 



27ri 



U e 



(4.8) 



with the integration path closed around z = 0. 

The function (t)-{z) is determined when ( |4.8| ) is substituted to Eq. ([4.6|). Note that 
the boundary condition A~(0, 77) = yl+(Lo, rf) is respected by Eq. 
Appendix A) is 



</)_(z) 



o2 Lo 



A 



The result (see 
(4.9) 



and 



(r? + Lo) 



(4.10) 



4.3 Analyzing the solutions 

Define 

Wo = \/4A^ 

The amplitude for 7* 7* scattering is obtained using Eq. ( p.25| ) 

4 TxT 



T 
and 



T{Qt Qi, s) 



T(Ql Ql s) 



lo { ujojln-^ In^ 



4 TxT 



1 



if 1^1 < Q^/s (4.11) 



^0 ( wojln^ In^ 



In- 



I^LWln^ln 



Mo 



if /^o > Q'A- (4.12) 



It is important to note that, although in the course of our derivation we seem to have 
lost the symmetry in Q| and Q2, the final result of the amplitude is fully symmetric again. 
In particular, the second term depends upon Q^/s only. The amplitude T~ reduces to T"*" 
when Q'^/s = Hq, i.e. when the dynamical infrared cutoff of the perturbative calculation 
reaches fi^, the limit of the nonperturbative infrared region. 

Let us consider, in some more detail, the s ^ 00 asymptotics for the case — ^ 
/Xq. We take s to be much larger than the Qf, but still within the region ( |2.10 ): 



1 < s/Q'^ < QVa^o- 



(4.13) 



In this region the asymptotics is obtained from the asymptotic behavior of the Bessel 
function /q 



T+(s ^ 00) 



4 ttt 



UJn 



27ruJo ln(s/Q2) 



s 



(4.14) 
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and the result is entirely perturbative. When s increases and eventually reaches the boarder 
line between and /~: 

we have to switch to T~ . With a further increase of s, initially, the second term in (|3.27| ) 
is not large and we can use the expansion of the Bessel function I2 for small arguments. 
In the asymptotic region 

QVmo <■ s/Q^ (4.16) 



the arguments of both Bessel functions are large. Expanding in both 



lns/Q2 



and in the 



ratio 



we find that in ( [4.17| ) the leading terms of both Bessel function terms cancel 



lns/Q2 ' 

against each other, and have to take into account first corrections. We obtain: 



r-(. 



00 



8 TxT 



a;2^(c^o ln(s/Q2))3/2 



I + ojoln^ + ^ \r? ^ + O 



^0 



In^ 
r)2 



(4.17) 



It is interesting to compare (4.14) and ( 4.17] ): the power behavior in s is the same 
in both regions. The difference lies in the preexponential factors: in the second region, 
we have a slightly stronger logarithmic suppression, and there is a logarithmic dependence 
upon the infrared scale /ig- 



4.4 DIS at low x 

Another case of interest is deep inelastic scattering on an almost real photon at very small 
X. This corresponds to the limit 



(4.18) 



and only the region / applies [Q = Qi): 
Tdis{Q^, s) = h[^^Q\ /in ^ In 



In^ V V /^o 



4 TxT 



(4.19) 



The Bjorken x is defined in a standard way: x = Q^/s. The flavor nonsinglet photon 
structure function is related to Tijjs via 



CJ- 



,7 7 
tot 



4 TT^ Ue 



dTpjsiQ^ s) 
din 



-^iVS'('''' Q 

71 - Ucm 71 - Ctery, S 4: 71 Uc^i U ill 

(4.20) 

Being aware of the fact that in the DIS limit (Q2 ~ A*o) saturation effects may 

be important [p4| we will not discuss them here at all. 
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We can consider two different asymptotic limits. The first one is In l/x ^ In Q^//iQ ^ 
1. In this hmit the structure function becomes: 



\xj OemW^ V2vra;o In'"' (1/x) V /^o . 



Eq. (4.21) gives the Regge hmit of the flavor nonsinglet structure function. Up to the 



preexponential factor this result agrees with the behavior of the flavor nonsinglet proton 



structure function found in [18|. 



Another asymptotic limit to be considered is 1 <C ln(l/x) <C ln(Q /^q) leading to 

Fls{x, Q") ~ X ^ ^ Vi-(V-) in(QV/^g) . (4.22) 

aemTTwg ^/2 ln(l/x) ln(Q2//^o) 



Eq. ( 4.22| ) comes from the asymptotic expansion of the first term in ( [4.19| ). The second 



term is subleading in this limit. Eq. ( [4.22 ) corresponds to the double logarithmic limit of 
the DGLAP equation. 

4.5 Scales 

It is instructive to compare our results for the high energy behavior of quark-antiquark 
exchange in 7* 7* scattering with those for gluon exchange, i.e. the LO BFKL Pomeron. 
For the latter it is well-known that, for sufficiently large photon virtualities and not too 
high energies, the internal transverse momenta are of the order of the photon virtualities 
and hence justify the use of perturbation theory (Fig. ^a). When energy grows, diffusion 
in In A;^ broadens the relevant region of internal transverse momenta, which has the shape 
of a "cigar" . Its mean size is of the order Vln s and eventually reaches the infrared cutoff 
/iQ. From now on, the BFKL amplitude - although infrared finite - becomes sensitive to 
infrared physics, and some modification due to nonperturbative physics has to be included. 




Figure 6: Energy dependence of the region of integration in a) BFKL and b) quark ladder. 



With the results of our analysis we now can make an analogous statement about quark- 
antiquark exchange. Since internal transverse momenta range between max{Q^ / s, /Xq) < 



k"^ < s or, equivalent ly 



max(- In s/Q^ In /xg/Q2) < In A;VQ^ < In s/Q^ 



(4.23) 



there is, for not too large energies, again, a limited region where transverse momenta stay 
above the infrared scale, and the use of perturbative QCD is justified. With increasing 
energy, the /c^-region expands and eventually hits the infrared cutoff /Xq. From now on 
the high energy behavior starts to depend upon infrared physics and requires suitable 
modifications. 

In order to understand, in the fermion case, the region of internal integration in more 
detail, let us first return to our amplitude, (illustrated in Fig.^). It satisfies the two 
dimensional wave equation ([4.2p. It is instructive to introduce the new variables 



t 



In- 



(5s 



which leads to the two-dimensional Klein-Gordon equation: 



i) 



1 



(4.24) 



Ut2 



4A 







(4.25) 




t 



Figure 7: Light cone variables for 
the wave equation. 



(here i\/4 A plays the role of the mass). 

As seen from ( [4.1[ ), our amplitude satisfies the bound- 
ary condition vl^(0, 0) = 1, and its normal derivatives 
vanish on the 'light cone' ^ = or 77 = 0. Let us 
illustrate this in Fig. |^: the upper end of the ladder 
we associate with the values /3 f« Q^/s and k^ Q^, 
which translates into ^ = 7; = 0ort = 2; = 0. 
At the lower end, when coupling to the external pho- 
ton, we put /? = 1, fc^ = Q^, which is equivalent to 
= r/ = In s/Q^ or t = In s/Q^, z = 0. Corre- 
spondingly, in Fig.|^ the t-axis points downwards: our 
ladder diagrams for 7* 7* describe the evolution (in t) from the initial point (tj, Zi) = (0, 0) 
(upper photon) to the final point {tf, zj) = (In s/Q'^, 0) (lower photon). As seen from 
( |4.1| ), all paths lie in the region of positive ^ and rj, i.e. they stay inside the square described 
by < 1^ < = In s/Q'^, < r] < r]f = In s/Q'^. Note the width in z: it starts from 
zero, increases up to its maximal value In s/Q'^, and finally it shrinks down to zero again. 

Let us confront this with the BFKL Pomeron: in Fig.^^ we have, once more, drawn 
the square which illustrates the internal region of integration. Apart from the difference 
in shape ("diamond" versus "cigar"), the most notable difference is the width in InA;^. 
In the fermion case it grows proportional to In s/Q^, i.e stronger than in the BFKL case 
where In k'^ grows as y^ln s/Q'^: the BFKL diffusion is replaced by a linear growth in the 
z-direction. Also the definition of 'internal rapidity' is different: in the BFKL case the 
vertical axis can be labeled simply by In 1//3, whereas in the fermion case our variable is 
t = In (/3 s/y^A;^ Q^) (in both cases, the total length grows proportional to Ins). 
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5. InfraRed Evolution Equation 



In this section we rederive our previous result from the IREE. We define the partial wave 
through the ansatz:: 

T{Ql Ql s) = l^^[4,y T(Q?//i^ Qi/A^ c.) . (5.1) 

For a more detailed description of the method we refer the reader to the original work 
1 17] as well as to some applications in Ref. The scale fl is introduced as 

an auxiliary infrared cutoff parameter, which, a priori, does not have to coincide with the 
infrared cutoff scale of nonperturbative physics, /xq. Before applying the IREE formalism 
to 7* 7* scattering, we have to discuss its applicability. Based upon our previous analysis, 
we have to distinguish between two cases. For any t-channel quark line we have the 
condition that /Iq < sa (3 = k"^ which, in the derivation of the IREE, will be replaced 
by /i^ < saf3 = k"^. On the other hand, we have the constraints that Q'l/s < f3, and 
Q2/S < a. Combining these conditions we have max{il'^, Q\ Q'^/s) < sa P = k'^, and this 
distinction (after identifying the auxiliary parameter /2 with the physical infrared cutoff ^q) 
leads to the the two regions, and . In the latter case, all internal transverse momenta 
are cutoff by /i, and this is the case where the IREE directly applies. In the former case, 
the internal momenta do not reach the value jl, i.e. the amplitude is independent of /2. 
This case, therefore, has to be studied separately. 

Let us begin with the genuine IREE region, 

Ql Ql/s < fi?. The IREE for the ampli- 
tude T is obtained by differentiating ( |5.1| ) with respect to In /i^ . In space this essentially 
corresponds to cutting of the diagram on two parts (see Fig P). 



F(yj) 



FCy,) 



Figure 8: The IRRE for 7* 7*. 

The equation reads 

+ — — + — — = 2NcFnc 3-^ (5.2) 

dyi dy2 Svr^ 
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with yi,2 = ln(Qf 2/A )> 



Fiuj, y) = /o(t^) e ^ ^ /or y > , (5.3) 



and Co = 47raeme(^) • 

The function /o was introduced in Ref. [17| 



with the property: 

/o(a; - ^) = 2vr2^2- (5-5) 

The function F{u}, y) is a sohition of the IREE for a quark-photon ampHtude |17, [l^]. The 
coefficient Co contains information about couphng of the photon to the ladder. Note the 
appearance in (5.2) of the factor 2NcFnc which incorporates a proper normahzation for 
the case of 7* 7* . 

We should impose boundary conditions for the case when only one photon is (nearly) 
on shell. For the boundary yi > 1/2 = it reads: 

(yi, y2 = 0) H = 2iVcF„c • (5.6) 

oyi 8 7r"= 

A similar condition is imposed at the second boundary y2 ^ yi = 0. When both photons 
are (nearly) on shell we have an additional equation 

^T-(yi = 2/2 = 0) = 2jV,F„ /'^^'^^= ^\ (5.7) 

We will find now a solution of Eq. ( |5.2D with the boundary conditions ( |5.6D and (p.7|). 
Note that the solution to be found is unique. 

A general solution of the inhomogeneous equation (|5.2| ) can be found as a sum of 
a particular solution (Tps) and a solution of the homogeneous equation (The), such that 
T~ = Tps + The- We will search for Tpg using the following ansatz, which is motivated by 
the right hand side of (|5.2|): 



fp,{uj) = f°,(w)e-('^--^o/*"')(^i + ^2) . (5.8) 
Substituting it to (U) we find ^^^(a;): 



2 



Jp.(u;) - «^_^^/4^2e (5.J) 
with 

_ Co(^)Co(^)2iV,Ke _ TTT _ TTT ^ 

8vr2g4c2 -54^2-4^4^4- l^-iU; 

Our goal now is to find The which satisfies the homogeneous version of Eq.( |5.2[) . Define 
two new variables: w = yi + y2 and u = yi — y2- Then 

dT 

cufheiio) + 2-^ = (5.11) 
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with the solution 

fUu^,u,w) = rO,(a;, n) = f^^ioj, - ys) (5.12) 

The boundary conditions (|5.6D and ( |5.7| ) will be used in order to find T^^. When T~ = 
Tps + The is substituted in Eq. ( ^.61 ) we get (?/2 = 0) 



2 ayi uj - /o/47r^ 



With the help of (p.5|) the rhs can be simplified: 



1^/?. + =- e-(-/"/-^)-/2. (5.14) 

2 5yi u - /o/47r^ 

Eq. (5.13) has a solution which, again, is given by a sum of a particular solution (.Rp^) and 
a solution of the homogeneous equation [Rhe), such that T^^ = iip^ + -R/ie- The particular 
solution is 

R^siu:, yi) = -^^4lt-2 (5.15) 
cj - /o/47r^ 

while the solution of the homogeneous equation is 

Rhe{u:,y{) = Rte-^y^'\ (5.16) 
Finally we use the boundary condition (5.7) in order to determine R^^,'- 

(5.17) 

UJ 

Combining all the results together we find the desired solution: 



+ 



{uj - /o/47r2) 



g-(y2-yi)/o/8 7r2 _ 

UJ 



which we can also write as 



(5.18) 



f-{uj) = -IGvr^K 



{uj - y/uj^ - a;g)2 _ /o/8 ,^2) (yi + y2) 



2 _ ,.,2 




16 TT^ KUJq 



, - u){yi + y2)/2 p (y2 - J/i) - "^0/2 p-J/i^ 



CJ' — UJn 



UJ 



(5.19) 



In ( 5.18 ) and ( 5.19 ), the first term on the rhs is the particular solution of Eq. (5^), whereas 
the second term is a solution of the homogeneous equation. Since we used the boundary 
condition ( |5.6| ), the solution ( 5.18|) (or ( ^.19[ )) is valid for yi > y2. Would we have started 



-20- 



with the boundary conditions at yi = 0, we would have ended up with ( |5.18| ) or ( 5.1S| ), 
but with yi and 2/2 being interchanged. 

At first sight, the distinction between the two cases yi < y2 and 7/2 < Hi might look 
somewhat strange. Returning with our result into putting /l^ = ^Ug, we reproduce 

exactly the expression ( [4.12 ) for the amplitude T~ . It is not difficult to identify (after 
integration over to) the first term on the rhs of (5.19) with the second term on the rhs of 
( 4.12D (the Bessel function I2), whereas the "nonsymmetric" second term equals Iq — 1. 
The final result is fully symmetric with respect to yi 2/2 interchange. The apparent 
puzzle is resolved because of the two-sheeted structure of the w plane for the partial wave 
representation of the Bessel function Iq (see Appendix B). 

Finally, we return to the other region, /i^ < Q\Q\js. As we have argued above, in 
this region the amplitude does not depend upon /I^. Starting from the ansatz (^), we 
have to require that the total derivative with respect to vanishes, i.e. 



a; 



r+ + 



dyi 



+ 



dy2 



0. 



(5.20) 



A solution to this equation has to be of the form ( 5.12| ), where the dependence upon the 
difference yi — y2 has to be fixed by boundary conditions. We require that, for /i^ = 
QIQ2/S the solution matches the solution ( ^.19 ), found for the other region, . Since, 
on the rhs of ( |5.19| ), the first term after integration over u> is proportional to the Bessel 

1, we are left with the 



function l2{uJo yln{sfl'^ /Q"^) ln{s/fl'^)) and vanishes for sfi^/Q"^ 
second term only and our result is 



-uj{yi+y2)/2 p(s/2-S/l) V'^^-'^o/2 



-yi ^ 



(5.21) 



After integration over w, this becomes /q — 1, in agreement with (4.11). 

In summary, we have been able to derive, within the IREE approach, both our solutions 
T~ and . This result allows for the possibility of studying also non-ladder diagrams 
which appear in odd signature exchange amplitudes. 



6. Reggeon Green's function 

A third method of describing the double logarithmic high energy behavior of quark- 
antiquark exchange has been developed in Ref. [^]. Kirschner has shown that the double 
logarithmic situation can be described in the same way as the BFKL Pomeron, namely by 
a two-reggeon Green's function with a kernel which is conformal invariant and holomor- 
phic separable. In particular, in the integral equation of the reggeon Green's function the 
transverse momentum integration has no explicit infrared or ultraviolet cutoff; therefore it 
has to be compared with our result for ^0 = 0. In the following we will show that, starting 
from Kirschner's Green's function, we arrive at the same result as in section 4. 



The reggeon Green function in Ref. |21| has been obtained by deriving and solving 



a linear BFKL- type equation for the quark exchange amplitude |20]. The Green function 
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has the form 



with 

Vt{u;, V, n) = 4V'(1) - i'i-iv + M + |) _ ^{iy + M + |) _ 

- iz. - ^ + -) - V(l + - y + 2) (6-2) 

and ^^{z) = ^ Inr(z). The function G is normahzed such that at the Born level it equals 
Go = -^jpr <^^(^± — ^±)) which is a propagator of two quark lines. 

At high energies the leading contribution comes from n = 0. After expansion of for 
small values of oo and u the Green functions reduces to 

1 1 + Av"^ fk'^Y" 

""('^''^'"^ = ^fcpiy-oo^' ^(^^-u.g + 4.2) (^J (6-3) 

The 7* 7* scattering amplitude can be obtained from G by integrating over momenta 
k in the low and upper blocks in Fig. |9|. Note that the propagators of the t-channel quarks 
(l/fc^) are already included in G. In Ref. |^] they were assumed to be purely transverse 
implying that 

aPs<^kj_ and af3s<^k'j_. (6.4) 

The a and /3 integrations pick up the poles of the s-channel quark propagators. The only 
trace of these integrations is in the ordering ( |3.16| ) in k'^/P (and in k'^/a). This ordering 
leads to the following limits for /? and a: 

OLmax = min{l , kj_ / Q^} Pmax = min{l ,k\ / Q\} . 

The scattering amplitude has the form 

TZ' ^' {Ql Ql s) = TTT I k^ fj^^ ^ / fj^^ J 

duj f a(3s Y nfh T \ ^(i k\ 

uj G[K±, k±, uj) . (0.0) 



2 TT i V I ^± I I ^± I 

The overall coefficient ttt as well as the factor uj in front of G account for the proper 
normalization and can be deduced using the Born approximation. Note that, following the 
derivation of G in |21|, we assume that any potential infrared and ultraviolet divergences 
stemming from the k± integrations are regularized by the iv dependent factor. This is why 
the oj integration should be done only after the k± integrations are performed. 
After the a and (3 integrations we obtain: 



max Pmax ^ \ 1 i \ i n n\ 

G[k±, k±, iO) . (6.6) 



\k±\\k±\ 
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It is worth to notice that the 7* 7* scattering amphtude can be brought to the form in 
which G is sandwiched between the photon impact factors: 



rpT 7 



NcFnc 



2-K 



2tt 



du! 



2TTiuj V I ^-L I I ^-L I , 

X kl, G{k^, k^, uj) ^{Ql, kl, io). 



(6.7) 



The photon impact factor has the form 

kl, LO) = Co 



kl 



e 1 



kl 



+ e 



kl 



1 



(6.8) 



with Co defined in the previous section. Integrating in ( |6.6[ ) over k±_ , k_\_ we get 
Sttt f dio f°° , 1 



rpTT 



vr 



2'Ki 



dv ■ 



The v integration leads to 



7 



4:TTT 



diO 
2TTi 



s 



q2 \ -V''^2_j^2/2 



1 fol 

Ql 



-uj/2 



Finally, after the oj integration we arrive at: 



rpT 7 

^G 



Attt 



UJn 



t^o 



Qi 



(6.10) 



(6.11) 



The expression ( |6.11 ) coincides exactly with the amplitude ( 4.11| ), with the restric- 
tion /iQ < Q^/s being removed. Formally we can say that the Green's function approach 
should be compared to our previous treatment (in particular, the linear equation in section 
4) with = being put to zero. A striking difference between the two treatments lies in 
the integration over transverse momenta. From the discussion in section 4 it follows that 
all transverse momenta inside the ladder lie in the region Q^/s < k"^ < s. In the inte- 



gral equation for the reggeon Green's function |21], on the other hand, the Regge energy 
factor {s/k±k'j_)'^ serves as a regulator, and the transverse momentum integration can be 
extended down to zero and up to infinity. 

An easy way to see how this works is the following. Return to our two-loop integral in 
section 3 and compare the lower cell, k2, with the k integration in ( |6.5| ). The integration 
region is illustrated in Fig. |^: interchanging the order of integration over P2, ^2 , we see that 
the ^2 integration goes from Q^/s to s and splits into two pieces. ForQ^/s < < Q^, the 
P2 integral extends from Q'^/s to k'^/Q'^, and for < < s we have fc^/'S < P2 < 1- 
If we now replace, in Fig. the upper cell by the reggeon Green's function, including 
the energy factor {f32S / y^QH^)'^ , we see that, as long as Re w > 0, we can extend the (32 
integral down to zero and enlarge the transverse momentum integral to the whole interval 
[0,00]. 
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Nevertheless, from the discussion in section 3 we 
know that the double logarithmic high energy behav- 
ior comes from the restricted region Q^/s < k2 < s 
(more precisely, from the shaded region in Fig. ^3). 
In the Green's function approach, this restriction of 
internal integration appears as soon as we replace 
the integrand of ( |6.1D by the n = 0, small-z/ and 
uj approximation, (|6.3[) . Before the making the ap- 
proximation, there is no restriction on paths con- 
necting, in Fig. ^3, the starting point at the top with 
the end point at the bottom. In particular, there 
will be paths outside the marked square. Turning 
to the leading high energy behavior and doing the 
n = 0, small-i^ and to approximations, it is easy to 



G(k, k,a)) 



Figure 9: 

function. 



7* 7* via reggeon Green 



verify that (O), when coupled to the upper photon, 

is identical to our A'^ function^: is satisfies the two-dimensional wave equation, and the 
discussion of section 4.5 applies. In this sense one might say that approximation which 
leads from ( |6.lD to ( |6.3| ) plays the same role as, in the BFKL case, the saddle point ap- 
proximation. 



7. Numerical estimates 



The final goal of this project should be a confrontation of the obtained results with the LEP 
data. We are not ready yet to produce such a comparison since we still miss a significant 
theoretical contribution arising from the flavor nonsinglet exchange. In this section we will 
present some first numerical estimates related to the resummation of the quark ladder. 
Using Eq. (2.15) the flavor nonsinglet contribution to (7^^^^ can be computed from 



the elastic amplitude ( 4.11| ) and ( 4.12| ). In our numerical estimates we will drop the flavor 
factor Fns- the missing flavor singlet piece will be estimated to have the same functional 
form as the nonsinglet piece, and ( |2.15| ) - with the factor Fns being replaced by ^ - is 
assumed to represent the sum of flavor singlet plus flavor nonsinglet. 

First let us demonstrate numerically the effect of the two kinematical regions appearing 
in ( [4.11| ) and ( 4.12| ). Fig. 10 shows the 7*7* cross section as a function of rapidity 
Y = Ins/Q^ for equal masses Qi = Q2 = Q^ = IGGeV^. Up to y ~ 6 this corresponds 
to the LEP data region. For the fixed value of as we use Qs(Q^) ~ 0.24. The three curves 
show the dependence of the cross section on the nonperturbative scale /io- The solid line 
shows the (unphysical) case /ig = (the region /+), the dotted line is /ig = 0.5GeV^, and 
the dashed line is //g = 1.2GeV^. The points where the different curves come together 
correspond to s ji^/Q'^ = 1. To the left of these points we have the hard domain where 
the perturbative QCD calculation is fully reliable and does not depend upon the infrared 



""^It is important to note that, in (6.3), the contours of integration of the u and ui integrals, are not 



independent of each other. For example, one first does the u integral (at fixed Reto > loq), then closes the 
remaining ui contour around the cut at — tJo < < tJo 
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b 
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Y= In (s / Q2) 

Figure 10: gg contribution to o'l^^ for various values of \x^. 

scale ;Ug. Note that for /Xg = 0.5 Gel/^, almost all LEP data are within the hard domain. 
In this region we should expect that the secondary reggeon contribution is described by 
pQCD, and one should not add a further nonperturbative reggeon. 

Fig. 11 compares the ladder resummation with the box diagram contribution^. A 
significant enhancement is observed. The enhancement grows at higher energies and reaches 
a factor of ten at y ~ 10. For comparison we also show the nonperturbative reggeon 



(dashed line) ~ s 



-0.45 



In Ref. P] this contribution was added to the box diagram in 



order to fit the data. We believe that within the hard domain our resummed ladder should 
replace the contribution of the phenomenological reggeon. This can be qualitatively seen 



from the Fig. 11. 



8. Conclusions 

In this paper we have considered quark-antiquark exchange in 7* 7* scattering. For the 
flavor nonsinglet channel we have obtained a closed expression for the cross section a'^ ^ , 

^Only the leading logarithmic contribution is taken for the box diagram. The results thus obtained are 
somewhat larger compared to the ones based on the exact expression for the box. 
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Figure 11: Various contributions to cr, 



valid within the double logarithmic accuracy of pQCD. The result depends on the four 
scales relevant for the problem s ^ > Q2 > fJ-Q- The infrared cutoff fiQ is introduced 
into the analysis as the momentum scale below which nonperturbative physics starts. 

One of the central observations of this work is the role of this infrared cutoff in the 
high energy behavior of the total cross section. For large photon virtualities , Q2, there 
is a region in energy where the internal transverse momenta lie in the hard region, and 
the result does depend upon the infrared cutoff. In this region the perturbative analysis is 
reliable. When energy increases beyond the value s = QI/a^O' region of integration 
starts to extend into the infrared region; as a first step of handling this infrared dependence 
we introduce a sharp cutoff, /ig < k'^, and the results starts to depend upon Hq. We find, 
however, that this dependence is weak (via a logarithmic prefactor only). 

There is a close analogy between this dependence in quark-antiquark exchange and 
the diffusion in the BFKL Pomeron. In both cases, the relevant region of integration grows 
with energy; in qq exchange this growth is faster than in the BFKL case. As far as the ^0 
dependence is concerned, there is also an analogy with the non-forwardness t of the BFKL 
pomeron. By preventing the quark transverse momenta to penetrate the nonperturbative 
domain, both fiQ (in qq scattering) and a nonzero t (in the BFKL case) trigger a change in 
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the high energy asymptotics of the amphtudes: the pre-exponential factor in the amphtude 
changes from l/Vln s to l/(Vln s)^. 

As a byproduct of our analysis we have derived the low-x asymptotics of the DIS flavor 
nonsinglet photon structure function. The power dependence on x is the same as for the 
proton structure function |18|. 

The resummed quark ladder serves as a model for a 'perturbative secondary reggeon'. 
It is remarkable that the resulting intercept ujq ~ 0.5 is very close to the one known from 
the high energy phenomenology. The large intercept is due to the fact that the leading 
contribution is double logarithmic and ujq ~ -y/aj. The square root dependence of ojo on 
as somewhat reduces the uncertainty due to unknown value of as- 

The study of the quark ladder has an obvious phenomenological motivation. The 
LEP data on 7* 7* are at energies at which the quark box (QPM) still gives a dominant 
contribution to the cross section. We have shown that the gluon radiation leads to a 
significant enhancement of the quark box and hence needs to be accounted for. The quark 
box contribution dies fast with energy and is correctly expected to be of no importance 
for 7* 7* scattering at a NLC. In contrast, the pQCD reggeon receives an enhancement of 
about a factor of ten compared to the quark box, and potentially can still give a noticeable 
correction to the dominant pomeron contribution. 

Our analysis is incomplete, since we have not yet calculated the flavor singlet quark- 
antiquark exchange. This step is more involved due to the admixture of t-channel gluons 
(which are in helicity states different from the BFKL Pomeron). The situation is similar 
to the flavor singlet contribution to the polarized structure function gi |]l5|: here it is 
the antisymmetric e-tensor which, in the gluon t-channel, projects on the nonleading he- 
licity configuration. The evaluation of this contribution to 7* 7* is in progress and will be 
reported in a further publication. 

In our work, the 7* 7* scattering via quark ladder exchange was investigated by three 
different methods. First we have derived and solved a Bethe-Salpeter-type equation for the 
scattering amplitude. We have discovered the two distinct kinematical regions described 
above: in the first region the resulting amplitude is purely perturbative and is not influenced 
by the nonperturbative scale /io, whereas in the second region, the amplitude is infrared 
sensitive. The second method which was applied is the nonlinear IREE. Using the IREE 
formalism known from the literature and adjusting it to 7*7*, we have reproduced the 
results obtained from the linear equation. A strong advantage of the IREE formalism lies 
in the fact that it allows to generalize to odd-signature exchange which includes nonladder 
Feynman diagrams. Finally, the method relying on the reggeon Green function was used. 
Using this method we successfully reproduced our result in the perturbative region [fiQ = 
0). We do not know whether this method can be applied to the flavor singlet exchange. On 
the other hand, the Green's function approach allows to generalize to the case of nonzero 
momentum transfer, and it also may allow to extend our results to single logarithmic 



accuracy. It is claimed in Ref. [21| that the Green function G (|6.lD is obtained within a 
single logarithmic accuracy, which is beyond the double log accuracy adopted in the present 
work. 
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A. Appendix 

In this Appendix we derive the result ( [4. 91 ) for the function (/>_. The function 4>-{^) can 
be found if ( |4.8| ) is substituted to Eq. ([4.61). Instead we will solve two simpler equations 
which follow from Eq. ( |4.6| ). This will be enough to uniquely determine 
The first equation is the boundary condition 

^-(0, r/) = ^+(Lo, 7?), (A.l) 

which can be rewritten as 

dz^^{z)e'^' = / _e''^ + ^^o/^ (A.2) 
0+ J0+ z 

Eq. ( [A.2| ) implies 

(t)^{z) = le^^»/^ + 4)-{z) (A.3) 
z 

with (/)_ (z) being a function regular at the origin. 

The second equation, which we are going to use, is obtained by differentiating Eq. 
(|4^ ) with respect to ^' . We get 

= Xj^^ dfjA"{C',fj) (A.4) 

which reduces to 

[ ^0_(^)e«(- + V-) = (A.5) 

and has to fulfilled for any ^. Eq. ( A.5 ) is solved by substituting (A^) and expanding (p^ 
in a power series in z 

4>-{z) = ^c„z". 

n=0 

Eq. ( [A.5| ) can be rewritten 

f ^e«- + ^« + ^«)/^ = - VC„ / dzz"-ie«(^ + ^/^). (A.6) 

Expanding on both sides of Eq. ( [A.5| ) the exponents in powers of 1/z and integrating over 
z we obtain 

t» ™ jm-k cm+k+l °° t2m-n 

^ k\im-k)\{m + l)\ m\{m-n)\' ' ' 

m=OA;=0 ^ ' ^ ' n=Om=n ^ ' 
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The coefficients in front of a given power should be equal for any p. This leads to 



m=(p-l)/2 



L 



2 m— p+1 




{p — m — 1)1 {2m — p + 1)1 (m + 1)! 



n=0 



\{P+n)/2 



{{p + n)/2)\{{p-n)/2)\ ' 



(Ai 



Introducing a new summation index m = 2m — p + 1 casts the left hand side of ( [A.8D 
into the same form as the right hand side. We finally find 



Cn 



1 L 



71-1 



A(n- 1)! 



and (p- (z) 



(A.9) 



B. Appendix 

In this appendix we compute the amplitude T~ by the inverse Mellin transform of the 
expression ( ^.19|) . Let us compute the amplitude T~ as a sum of two separate terms 
= Ti + T2 where 



T^ = K 



duo 



Cc 



2ni \fx 



r,2 



UJ 



-4 ^ 1^- /o/Svr^ 



(B.l) 



T2 = 16 vr^ ujq h 



/i2 



q2\ a/^^~"o/2 



(u;-/o/47r2) V /" 



2 



1 



(B.2) 



As we will see below, the w-integration path Ccut goes around the square root branch cut 
from — (jQ to LiJo- 

For the first term in the amplitude we have 



"° (W - v^U;2 _ ^2)2 



- 



Denote by ^0 = ln(s /i^/Q'*) and r/o = ln(s//i2) and introduce a new complex variable 



(B.3) 



(B.4) 



The u; integral then turns into a contour integral in the z-plane, and the integration path 
encircles the origin z = 0. We obtain: 



Ti = KlGvr'^wg / 



dz z 
2m 



,rio wo z/2 + go 1^0/2 z 



KlGvr^ cj^ — l2(u;o \/Co%), (B.5) 



quite in agreement with the second term on the rhs of ( [4.12| ). 
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For the second part T2 we have 



■J' 



01 
Ql 



(B.6) 



Denote by ^1 = ln(s/Qf) and rji 



ln(s/(52)) introduce z as above. We obtain: 



To 







dz 



0+ 



2 vr i z 



(B.7) 

which agrees with the first term on the rhs of ( 4.12 ). 

Note that the rhs of (|B.6|) , although it looks asymmetric in yi and 2/2 > is completely 
symmetric under the the exchange yi ^ y2- This is because, in ( [B.?]) , we have the freedom 
to substitute z 1/z, i.e. to interchange the two momentum scales and Q^- The 
transformation ( [B.4D maps a single sheet z-plane to a two-sheeted uj plane. The simple 
change of variable: z 



1/z is equivalent to replacing lo — y'i 



UJn 



UJ + ^/up■ 



(and to moving the uj contour from the first to the second sheet). Also, the last term on 
the rhs of (|B.6|) leaves the freedom to multiply with (■^)^ or with (tIj)'^: in both cases, 
the oj integration gives 1. 



References 

[1] E. A. Kuraev, L. N. Lipatov, and F. S. Fadin, Sov. Phys. JETP 45 (1977) 199; Ya. Ya. 
Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 22. 

[2] P. Achard et al. [L3 Collaboration], Phys. Lett. B 531 (2002) 39. 

[3] G. Abbiendi et al. [OPAL Collaboration], Eur. Phys. J. C 24 (2002) 17. 

[4] S. J. Brodsky, F. Hautmann and D. E. Soper, Phys. Rev. D 56 (1997) 6957. 

[5] J. Bartels, A. De Roeck and H. Lotter, Phys. Lett. B 389 (1996) 742. 

[6] A. Donnachie, H. G. Dosch and M. Rueter, Eur. Phys. J. C 13 (2000) 141; Phys. Rev. D 
59 (1999) 074011. 

[7] J. Kwiecinski and L. Motyka, Eur. Phys. J. C 18 (2000) 343. 

[8] N. N. Nikolaev, J. Spcth and V. R. ZoUer, J. Exp. Theor. Phys. 93 (2001) 957, [Zh. Eksp. 
Tear. Fiz. 93 (2001) 1104]. 

[9] S. J. Brodsky, V. S. Fadin, V. T. Kim, L. N. Lipatov and G. B. Pivovarov, JETP Lett. 76 
(2002) 249, [Pisma Zh. Eksp. Tear. Fiz. 76 (2002) 306]; hep-ph/0111390. 

[10] M. Cacciari, V. Del Duca, S. Frixione and Z. Trocsanyi, JHEP 0102 (2001) 029. 

[11] V. S. Fadin and L. N. Lipatov, Phys. Lett. B 429, (1998) 127; 
G. Camici and M. Ciafaloni, Phys. Lett. B 430 (1998) 349. 

[12] J. Bartels, S. Gieseke, C.F. Qiao, Phys. Rev. D 63 (2001) 056014; Erratum- ihid D 65 
(2002) 079902; J. Bartels, S. Gieseke, A. Kyrieleis, Phys. Rev. D 65 (2002) 014006; J. 
Bartels, D. Colferai, S. Gieseke, A. Kyrieleis, Phys. Rev. D 66 (2002) 094017. 



-30- 



[13] V. S. Fadin, A. D. Martin, Phys. Rev. D 60 (1999) 114008; V. S. Fadin, D. Yu. Ivanov, M. I. 
Kotsky, Nud. Phys. B 658 (2003) 156. 

[14] M.G. Ryskin and A.G. Shuvaev, Eur. Phys. J. C 25 (2002) 245. 

[15] J. Bartels, B. I. Ermolacv and M. G. Ryskin, Z. Phys. C 72 (1996) 627. 

[16] V. G. Gorshkov, V. N. Gribov, L. N. Lipatov and G. V. Frolov, 5*01;. J. Nud. Phys. 6, 95 
(1968), [Yad. Fiz. 6, 129 (1967)]; Sov. J. Nud. Phys. 6 (1968) 262, [Yad. Fiz. 6 (1967) 
361]. 

[17] R. Kirschner and L. N. Lipatov, Nud. Phys. B 213 (1983) 122. Sov. Phys. JETP 56 
(1982) 266, [Zh. Eksp. Teor. Fiz. 83 (1982) 488]. 

[18] B. I. Ermolaev, S. I. Manaenkov and M. G. Ryskin, Z. Phys. C 69 (1996) 259. 

[19] J. Bartels, B. I. Ermolaev and M. G. Ryskin, Z. Phys. C 70 (1996) 273. 

[20] J. Kwiecinski, Phys. Rev. D 26 (1982) 3293; 
R. Kirschner, Z. Phys. C 31 (1986) 135. 

[21] R. Kirschner, Z. Phys. C 67 (1995) 459. 

[22] V. M. Budnev, I. F. Ginzburg, G. V. Meledin and V. G. Serbo, Phys. Rept. 15 (1974) 181. 
[23] A. V. Kotikov and L. N. Lipatov, hep-pli/0208220. 

[24] K. Itakura, Y. V. Kovchegov, L. McLerran and D. Teaney, liep-ph/0305332. 

[25] R. Kirschner, L. Mankiewicz, A. Schafer and L. Szymanowski, Z. Phys. C 74 (1997) 501. 



- 31 - 



